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We study the effect of uniaxial strain on the transmission and the conductivity across 
a strain-induced barrier in graphene. At variance with conventional studies, which consider 
sharp barriers, we consider a more realistic, smooth barrier, characterized by a nonuniform, 
continuous strain profile. Our results are instrumental towards a better understanding of the 
transport properties in corrugated graphene. 
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Graphene is a single layer of sp 2 carbon atoms, arranged as an honeycomb lat- 
tice. Its fabrication in the laboratory [l| immediately stimulated the interest of 
both the experimental and theoretical communities, and many applications, which 
initially could only be speculated, now appear feasibile. In particular, electronic 
quasiparticles in graphene are characterized by a band structure consisting of two 
bands, touching at the Fermi level in a linear, cone-like fashion at the so-called 
Dirac points ±K, and a linearly vanishing density of states (DOS) at the Fermi 
level @, 0] • This implies in this novel condensed matter system the possibility of 
Klein tunneling across barriers 0-0], i.e. perfect transmission across energy barri- 
ers, which was predicted in the context of quantum electrodynamics at relatively 
much larger energies. 

Graphene, like most carbon compounds, is also characterized by quite remark- 
able mechanical properties. Despite its reduced dimensionality, graphene possesses 
a relatively large tensile strength and stiffness with graphe ne sheets being 
capable to sustain elastic deformations as large as « 20% [llMl5l|. Larger strains 
would then induce a semimetal-to-semiconductor transition, with the opening of 
an energy gap lil-22T|. 

Recently, it has been suggested that graphene-based electronic devices might be 
designed by suitably tailoring the electronic structure of a^ graphene sheet under 
applied strain (the so-called 'origami' nanoelectronics) (23| . Indeed, a consider- 
able amount of work has been devoted to the study of the transport properties in 
graphene across strain- induced single and multiple barriers 241 ] . It has also been 
suggested that strain may induce relatively high pseudo- magnetic fields [25] , whose 
effects have actually been confirmed experimentally in graphene nanobubbles grown 
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on top of a platinum surface [26]. Indeed, the effect of the strain-induced displace- 
ment of the Dirac points in reciprocal space can be (formally) described in terms 
of the coupling to a gauge field. However, applied strain also induces a variation 
of the Fermi velocity, v-p. In particular, uniaxial strain implies a Fermi velocity 
anisotropy, while an inhomogeneous strain implies a nonuniform (i.e. coordinate 
dependent) velocity profile. This can also be realized in the presence of smooth 
potential barriers, where it has been demonstrated that a nonuniform space vari- 
ation of the underlying gate potential would result in a modulation of the Fermi 
velocity [27M29I] . Moreover, there is considerable evidence, both experimental [13] 
and theoretical [23], that barrier edge effects are also important to determine the 
transport properties across corrugated graphene. 

Close to the Fermi energy and in the unstrained case, the electrons dynamics is 
governed by the linearized Hamiltonian 

H = hvpcr ■ p, (1) 

where cr = (a x ,a y ) is a vector of Pauli matrices, associated with the in-plane 
spinorial nature of the quasiparticles in graphene. Eq. (prj can also take into ac- 
count intervalley processes K O — K, which however can be safely neglected, at 
sufficiently low energies. 



Applied strain is then described by means of the strain tensor 17] e = ^e[(l — 
u)i + {l + v)A{9)], where A{9) = a z e 2t0r7v . Here, e is the strain modulus, 9 the angle 



along which strain is applied, and v = 0.14 is the Poisson ratio for graphene [2d. I22l. 
EjlL I32I ] . Starting from a more general, tight-binding Hamiltonian [3] , and expanding 
to first order in the strain modulus, one obtains an anisotropic dependence on the 
strain ang le 9, already at linear order in the impulses & This can mapped back to 
a linear Hamiltonian as in Eq. ([T]) [22! ). where now impulses are reckoned from the 
shifted Dirac points, and the Fermi velocity is anisotropic and possibly coordinate- 
dependent [32j. 

We therefore consider a smooth strain barrier, characterized by a nonuniform, 
continuous strain profile e = e(£), with 

£ ® = tanh( J D/4a) (l + e"«/ a ~ l + e -«-£»/a) ' (2) 

where £ is the coordinate along the strain direction, forming an angle 9 with the 
crystallographic x axis. Such a strain profile is essentially flat for |£ — D/2\ <C a, 
where e(£) pa e 0; and for |£ — D/2\ ^> a, where e(£) « 0. In the limit a/D — > 0, 
Eq. ([2]) tends to a sharp barrier. The linear extent a, over which the strain profile 
Eq. ([2]) varies appreciably, is naturally to be compared with the lattice step 0, at 
the microscopic level, and with the Fermi wavelength Af = hvp / {2ttE), where E 
is the energy of the incoming electron. While a smooth profile can be expected on 
quite general grounds, the approximation of a sharp barrier is expected to hold 
well whenever a<o< Af, i-e. at sufficiently large incident energies. On the other 
hand, the details of the strain profile come into their own when a ~ X-p. 

Single electron tunneling, and thus the majority of the transport properties of 
interest, can then be inferred by solving the stationary Dirac equation associated 
to Eq. ([T]), now including the nonuniform strain, Eq. ([2]) [32]. For |£| — > 00, the 
solutions for the scattering problem are therefore known analytically [32I ] (and 
refs. therein). Integrating the scattering equations from large positive £ backwards 
to large negative £, and comparing with the known analytical solution, one may 
extract the reflection coefficient r, relative to an incident wave with unit amplitude 
incoming from £ > 0, as the Fourier weight with respect to its negative frequency 
component, whence the transmission T(E, <p) at given incidence energy E follows 
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Figure 1. (Color online) Left: Tunneling transmission vs incidence angle ip across a smooth strain barrier, 
with D = 100 nm, and incidence energy E = 167 meV (Ap = 0.6 nm). Right: Normalized conductivity, 
o a (E), vs incidence energy E across a smooth strain barrier, Eq. Ht . In both panels, eo = 0.01, 9 = tt/2 
(strain is applied in the armchair direction), while different curves refer to (a) sharp barrier, a = 0; (b) 
a = 10~ 2 D; (c) a = 10 _1 D. 

straightforwardly. 

Fig. Q] (left panel) shows the transmission T(E,ip) as a function of the incidence 
angle ip across strain-induced sharp and smooth barriers, Eq. ([2]), with strain ap- 
plied along the armchair direction (0 = tt/2). One observes that, upon increasing 
the smoothing parameter a/D, the oscillations, characteristic of Klein tunneling 
across energy barriers in graphene, get damped, while their envelope (lower bound) 
increases. The dependence of the transmission T(E, ip) on the incidence angle ip is 
only apparently asymmetric, as we are restricting to quasiparticles with momen- 
tum centred around a given Dirac cone, say +K. Symmetry T(E, (p) = T(E, —ip) 
would be restored when the effect from the neighbourhood of both Dirac cones is 
included. 

The conductivity can then be str aigh tforwardly related to the transmission by 
means of the Landauer formula [33|, |34| , as 



a(E) = a D 



E 



1 (E,(p) cos ip 

-tt/2 



2vr' 



(3) 



where <7o = 4e 2 //i is twice the conductance quantum, and the conserved component 
of transmitted momentum, i.e. that parallel to the barrier, has been related to 
the incidence angle through k y = E/{hvp)sva.ip. In Eq. ([3]), only the propagating 
modes have been included in the integration. One is then prompted to define the 
adimensional conductivity 



& a (E) 



°{E) 



(4) 



Fig. [T] (right panel) shows the reduced conductivity, Eq. (|3J), as a function of 
incident energy E, for tunneling across sharp and smooth barriers, Eq. ([2|). One ob- 
serves Fabry-Perot oscillations, whose amplitude is reduced by increasing smooth- 
ing (i.e., increasing a), the barrier then tending to be a more regular function. 
Also, the overall increase of the transmission is reflected in an enhancement of the 
conductivity. 

In conclusion, we have studied the effect of nonuniform strain on the conduc- 
tivity across smooth strain barriers in graphene. While an increase of smoothing 
reduces the oscillations of the transmission as a function of the incidence angle, 
one finds a reduction of the Fabry-Perot oscillations and an overall enhancement 
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of the conductivity as a function of incidence energy. These results should help 
understanding the properties of corrugated graphene. 
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